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ABSTRACT 

The lensing cross section of triaxial halos depends on the relative orientation between a halo's principal axes 
and its line of sight. Consequently, a lensing subsample of randomly oriented halos is not, in general, randomly 
oriented. Using an isothermal mass model for the lensing galaxies and their host halos, we show that the lensing 
subsample of halos that produces doubles is preferentially aligned along the lines of sight, whereas halos that 
produce quads tend to be projected along their middle axes. These preferred orientations result in different 
projected ellipticity distributions for quad, doubles, and random galaxies. We show that w 300 lens systems 
must be discovered to detect this effect at the 95% confidence level. We also investigate the importance of halo 
shape for predicting the quad-to-double ratio and find that the latter depends quite sensitively on the distribution 
of the short-to-long axis ratio, but is otherwise nearly independent of halo shape. Finally, we estimate the 
impact of the preferred orientation of lensing galaxies on their projected substructure mass fraction, and find 
that the observed alignment between the substructure distribution and the mass distribution of halos result in a 
negligible bias. 

Subject headings: galaxies, halos, lensing 



1. INTRODUCTION 



Statistics of lensing galaxies have been used as cos- 
mological and galaxy formation probes sin ce early in 
the m odern history of gravitational lensing dTurner et alj 
Il984l) ^ Lensing rates can be used to constrain dark en- 



ergy (Pukugita et aDfT99l IChad 120031: iMitchell et al.ll200l 
ICha&.2007; ,Og uri et al.ll2007h. to probe the struc ture of lens- 
inggalaxies jKeetonl 120011; iKochanek & White! 12001; Chai 



2005b. and to probe galaxy evolution ("Ch ae & Mag i2003l ; 



Ofek et al.| [2003; Rusin & Kochanek 2005). While the use of 



lensing statistics as a cosmological probe has had mixed suc- 
cess, particularly early on, it remains a unique probe with en- 
tirely different systematics from more traditional approaches. 
Consequently, lensing statistics are likely to remain a funda- 
mental cross-check of our understanding of cosmology and 
galaxy evolution. 

One of the difficulties that confronts the study of lensing 
statistics is that, in general, the halo population that produces 
gravitational lenses can in fact be a highly biased subsample 
of the general halo population. For instance, it has long been 
known that while early type galaxies compose only « 30% 
of all luminous galaxies, the majority of lensing galaxies are 
in fact early type since these tend to be more massive and re- 
side in more massive halos than their late counterparts. By the 
same token, lensing early type galaxies tend to have higher lu- 
minosity and velocity di spersions than non-lensing early type 
galaxies dMoeller et alj [2006 ; Bolton et alj l200d) . Overall, 
then, when interpreting lensing statistics, one ought to always 
remember that by selecting lensing galaxies one is automati- 
cally introducing an important selection effect that can signif- 
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icantly bias the distribution of any galaxy observable that has 
an impact on the lensing probabilities. Here, we consider one 
such source of bias, the triaxiality of galaxy halos.* 

That halo triaxiality can have important consequences for 
lensin g statistics has been kno wn for several years. For in- 
stance JOgurT&^^eto^ (|2003) have shown that triaxiality can 
significantly enhance the optical depth of large image sep- 
aration lenses. Similar conclusions have been reached con- 
cer ning the formation of giant arcs by lensing clusters (see 
e.g. lOguri et alj2003l;lR"ozo et alj200"6allHennawi et al.l2007l 
and references therein). Curiously, however, little effort has 
gone into investigating how observational properties of lens- 
ing galaxies can be different from those of the galaxy popula- 
tion as a whole due to the triaxial structure of galactic halos. 
This work addresses this omission. 

The first observable we consider is the projected axis ratio 
of lensing galaxies. Roughly speaking, given that non-zero 
ellipticities are needed in order to produce quad systems, one 
would generically expect lenses that lead to this image con- 
figuration to be more elliptical than the overall galaxy popula- 
tion. Likewise, lensing galaxies that produce doubles should, 
on average, be slightly more circular than a random galaxy. 
There can, however, be complications for these simple pre- 
dictions due to halo triaxiality. For instance, given a prolate 
halo, projections along the long axis of the lens will result 
in highly concentrated, very circular profiles. Will the in- 
crease in Einstein radius of such projections compensate for 
the lower ellipticity of the system, implying most quads will 
be projected along their long axis, or will it be the other way 

* Throughout this work, we will be using the term galaxy and halo more or 
less interchangeably. The reason for this is that we are primarily focused on 
the impact of halo triaxiality on the lensing cross section, and the latter de- 
pends only on the total matter density. Consequently, differentiating between 
halo and galaxy would only obfuscate presentation and introduce unneces- 
sary difficulties. For instance, while modeling the total matter distribution as 
isothermal is a reasonable approximation, neither the baryons nor the dark 
matter by itself is isothermally distributed. Thus, it is much simpler to adopt 
an isothermal model, and refer to the baryons plus dark matter as a single en- 
tity, than to try to differentiate between the two. Likewise, when discussing 
triaxiality, what is important in this work is the triaxiality of the total matter 
distribution. 
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around? Clearly, the relation between ellipticity and lensing 
cross sections is not straightforward once triaxiality of the 
lensing galaxies is taken into account, but it seems clear that 
there should be some observable difference between the el- 
lipticity distribution of lensing galaxies and that of all early 
types. Interestingly, no such difference ha s been observed 
dKeeton et alJll997l: lRusin & Tegmarkl2001h . which seems to 
fly in t he face of our expec tations (though see also the discus- 
sion in lKeeton et alJfT998h . Is this actually a problem, or will 
a quantitative analysis show that the consistency of the two 
distributions is to be expected? Here, we explicitly resolve 
this question, and demonstrate that current lens samples are 
much too small to detect the expected differences. 

Having considered the ellipticity distribution of random and 
lensing galaxies, it is then a natural step to investigate the im- 
pact of halo triaxiality on predictions of the quad-to-double 
ratio. Specifically, it is well known that the quad-to-double 
rati o is sensitive to the e llipticity distribution of lensing galax- 
ies (iKeeton et alJll997l) . so if lensing can bias the distribu- 
tion of ellipticities in lensing galaxies, then it should also 
affect the predicted quad-to-double ratios. This is an im- 
portant point because it has been argued that current predic- 
tions for the quad-to-double ratio are at odds with observa- 
tions. More specifically, the predicted quad-to- double ratio 
for the CLASS (Cosmic Lens All-Sky Survey, iMvers et al] 
l2003t iBrowne et al. 2003) sample of gravitational lenses is 
too low relative t o observations (iRusin & Tegmarlj 120011: 
iHuterer et alj|2005h . Curiously, however, recent work on the 
quad-to-double ratio observed in the SQLS (Sloan Digital Sky 
Survey Quasar Lens Search, lOguri et aP l2006t llnada et aH 
12007). suggests that the exact opposite is true for the latter 
sample, namely, theoretical expectations are too high relative 
to observations (IOgurill2007l) . In either case, it is of interest to 
determine how exactly does triaxiality affects theoretical pre- 
dictions, especially since the aforementioned difficulties with 
the CLASS sample has led various authors to offer possibil- 
ities as to how one might boost the expected quad-to-double 
ratios. Specifically, one can boost the quad-to-double ration 
in the class sample either from the effect of massive satellit e 
galaxies near the lensing galaxies (ICohn & Kochane5l2004l) . 
or through the large-scale e nvironment of the lensing galaxy 
(IKeeton & ZabludofjIIOOl . Clearly, we should determine 
whether halo triaxiality can be added to this list. 

This brings us then to the final problem we consider here, 
namely whether the substructure population of lensing galax- 
ies is different from that of non-lensing galaxies. Specifi- 
cally, we have argued that lensing galaxies will not be isotrop- 
ically distributed in space. Since the substructure distri- 
bution of a dark matte r halo is typically aUgned with its 
parent halo 's long ax is (Zentner et al."2005'; 'Libeskin d et alJ 
llOOS; Agus tsson & B rainerd 2006; Azzaro et al. 2006), it fol- 
lows that the projected distribution of substructures for lens- 
ing galaxies may in fact be different for lensing halos than 
for non-lensing halos. Such an effect could be quite im- 
portant given the claimed tension between the Cold Dark 
Matter (CDM) p redictions for th e substructure mass frac- 
tion of halos (se e Mao et al. 2004) and their observe d values 
(iDalal & Kochanek 2002a : Kochanek & Dalalll2004l) . Like- 
wise, such a bias would impact the predictions for the level 
of astrometric and flux perturbations pro duced by dark mat- 
ter substructures in gravitational lenses (iRozo et al.l l2006bl: 
IChen et al.ll2007h . Here, we wish to estimate the level at 
which the projected substructure mass fraction of lensing ha- 
los could be affected due to lensing biasing. 



The paper is organized as follows: in section |2] we derive 
the basic equations needed to compute how observable quan- 
tities will be biased in lensing galaxy samples due to halo tri- 
axiality. Section [3] presents the model used in this work to 
quantitatively estimate the level of these biases, and discusses 
how lensing halos are oriented relative to the line of sight as 
a function of the halos' axes ratios. Section|4]investigates the 
projected axis ratio distributions of lensing versus non-lensing 
galaxies, and demonstrates that present day lensing samples 
are too small to detect the triaxiality induced biases we have 
predicted. Section|5]discusses the problem of the quad to dou- 
ble ratio, and section |6] demonstrates that halo triaxiality bi- 
ases the projected substructure mass fraction in lensing halos 
by a negligible amount. Section [T] discusses a few of the ef- 
fects we have ignored in our work and how these may alter 
our results, and finally section [8] summarizes our work and 
presents our conclusions. 

2. LENS BIASES INDUCED BY TRIAXIALITY 

We begin by deriving the basic expressions on which we 
rely to estimate the effects of halo triaxiality on the observed 
properties of lensing galaxies. In particular, we show that 
since the lensing cross section for triaxial lenses is in gen- 
eral not spherically symmetric, this implies that a population 
of randomly oriented halos produces a non-random lens pop- 
ulation. Finally, we show that the induced non-randomness of 
the lensing halo population can alter the mean observational 
properties of these halos relative to the general halo popula- 
tion. 

2.1. The Lensing Cross Section 

Let p be a set of parameters that characterizes the projected 
gravitational potential of a halo. For instance, p can be the 
Einstein radius of the lens, its ellipticity, and so on. Given a 
background source density n^izs) and a halo density nh{p,Zh), 
and in the absence of a flux limit, the mean number of lensing 
events per unit redshift per area is given by 

, , ,c^ = n.sfc)nA(Pi^'i):i— cr(p,Z/,,Z.O (1) 
dzsdzhdil dzsdzh 

where x is the comoving distance to the appropriate halo or 

source redshift, and 

CT(p,2/i,Zi)= / '^"y- (2) 

J lensing 

The integral is over all regions of the source plane that pro- 
duce lensed images of interest. For instance, if one were in- 
terested in quadruply imaged sources, the integral would be 
over all source positions that result in four image lenses. The 
quantity a is called the lensing cross section, and of particu- 
lar interest to us will be the cross sections cr^'*'' for producing 
A^-image systems. 

In reality, one always has some flux limit F,„i„ which cor- 
responds to a minimum source luminosity L„,i„. Fortunately, 
the above argument is easily generalized: let dns{L,Zs)/dL be 
the number density of background sources with luminosity L. 
Then, the mean number of lensing events becomes 



dNie 



dx dx I ,2 
■nh- ;— / d y 



dn,{L,Zs) 

dL — . (3) 

dL 



dzsdzhdfl dzs dzh j JL„„„/fi(y) 
If the source luminosity function can be approximated by a 
power law dns{L,Zs)/ dL = AL'" (note both A and a can de- 
pend on Zs), the above expression reduces to 



dNi 



lens, 



dx dx 



J J ,n,=nb{> L„n„)ni, — —aB(p,a,Zh,Zs) (4) 
azsdzhdu dZsdzh 
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where nh(> L,ni„,Zs) is the number density of sources above 
the flux limit in the absence of lensing, and as is given by 



crB(p,ZA,z.5,a)= / d y pi(y) 



\a-\ 



(5) 



where iijy ) is the total magnific ation of a source at position y. 
Following iHuterer et alj ( l2005h . we call as the biased cross 
section. Indeed, since the distribution of magnifications pip) 
among all lensing events is given by 

/'(a^) = ~ j "^"y '5(M(y) - (6) 

where a is the (unbiased) lensing cross section defined in Eq. 
121 then we can rewrite Eq. |4]as 



where 



(m"-') = j dfi p(m)m 



a-l 



(7) 



(8) 



Thus, the net effect of gravitational magnification on the fre- 
quency of lensing events can be summarized as a biasing fac- 
tor (/^"~') that multiplies the unbiased lensing cross section 
a. 

2.2. Triaxiality and Lensing Biasing 

Let P characterize the mass distribution of a triaxial halo, 
and let n be the orientation of the halo's long axis relative 
to the line of sight. The halo's two dimensional potential is 
then characterized by a new set of parameters p(P, n) which 
depend on the halo properties P and the particular line of sight 
n along which the halo is being viewed. For instance, the 
vector P can include such halo properties as halo mass and 
axis ratios, whereas p could include parameters such as the 
Einstein radius of the projected mass distribution as well as 
the projected axis ratio. 

As discussed above, the mean number of lensing events per 
unit redshift by a halo along a given line of sight is given by 
Eq. |4] For convenience, we define the halo and source surface 
densities dHh/dP and dY^ Jdzs via 



dPdzh 



dnh dx 
dP dzh 



d'^s _ dx 

dZs dZs 



(9) 



(10) 



In terms of these surface densities, and assuming a randomly- 
oriented distribution of halos, the mean number of lenses per 
unit area as a function of their orientation n is given by 



dNu. 



1 c/Ej dY,i, 



dPdndZsdzhdfl 2tt dZs dPdzh 



crB(p(P,n),z/,,z,,a). (11) 



The prefactor of l/(27r) arises from the fact that dn^/dPdn = 
{dnij dP) / {111) due to our assumption of randomly oriented 
halos. ^ We emphasize that Eq. [TT] characterizes the number 
of lenses as a function of the relative orientation n between 
the halo's major axis and the line of sight. Thus, to compute 
the total number of lenses irrespective of halo orientation, we 

' If n denotes the angle between the hne of sight and a specified halo 
axis, and given that n and — n correspond to the same line of sight, then it is 
evident that the space of all fines of sigfit is simply 5^/Z2 - a sphere with its 
diametrically opposed points identified. The volume of such a space with the 
usual metric is thus simply 2it. 



would simply integrate the above expression over all lines of 
sight n. 

There is an absolutely key point to be made concerning Eq. 
nn which provides the motivation behind this work. Specif- 
ically, we note that the number of lenses is proportional to 
o'B(p(P,n)). This implies that even though the overall halo 
population does not have a preferred orientation in space, the 
lens population is not randomly oriented, a fact which can 
have observable consequences. In particular, given an observ- 
able halo property /(P,n) that depends on the line of sight 
projection (e.g. the projected axis ratio or projected substruc- 
ture mass fraction), the mean value of / over all P halos is 
simply 



d-n 



/(P,fi), 



(12) 



whereas the mean value of / over all lenses is given by 
1 



(/|p; 



lenses 



d~h 



I IP\ / ^'^B(P(P,fl))/(P,fl) (13) 

(<Tb P) ./ l-K 



where (ctb) is the average value of as over all lines of sight, 

(14) 



{^b\P)= / ^aB(p(P,n)). 

Thus, in general, one expects that the mean value of / over 
all lenses and over all halos will be different. In the next few 
sections, we identify a few halo properties that depend on line 
of sight projection, and determine whether lensing biases in- 
duced by triaxiality are likely to be significant. 

3. THE MODEL 

We estimate the impact of halo triaxiality on the proper- 
ties of lenses by considering a triaxial isothermal profile. The 
merit of this approach is its simplicity: because of the simple 
form of the matter distribution in this model, we can compute 
all of the relevant quantities in a semi-analytic fashion, and the 
main features of the model can be easily understood, thereby 
providing an important reference point for investigating more 
elaborate models. Moreover, by working out in detail a sim- 
ple analytic model, our results provide an ideal test bed for 
more involved numerical codes, which would then allow us to 
investigate how our conclusions are changed as more compli- 
cated models are allowed (Chen et al. 2007, in preparation). 

3.1. Semi- Analytical Modeling 

Our analytical halo model is that of a simple triaxial isother- 
mal profile of the form 



p{x) = N{quq2) 



1 



InGx^/qj + y^ + Z^/qj 



(15) 



where qi and q2 are the axis ratios of the profile and we have 
chosen a coordinate system that is aligned with the halo's 
principal axes, and such that 1 > > qi-'^ The normalization 
constant A^(^i , qi) is chosen to ensure that the mass contained 
within a sphere of radius r be independent of the axis ratios 
for fixed velocity dispersion <t^, the latter being the velocity 
dispersion of the Singular Isothermal Sphere (SIS) obtained 
when qi =q2 = I. 

Let then (6',0) denote a line of sight. In appendix iBl we 
show that the corresponding projected surface mass density 

i.e. qi is the ratio of medium to long axis of the halo, whereas q2 is the 
ratio of the short to long axis. The motivation behind our particular choice of 
axis labeling will be made clear momentarily. 
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profile is that of a Singular Isoth ermal Ellipsoid (SIE) which, 
following lKormann et al.l (Il994h . we write as 



S.2 



2G 



(16) 



where both q and a,, are known functions of , q2 and, in the 
case of (T,,, of A^(^i , qojf^y (see AppendixlElfor details). In the 
above expression, a and q are the effective velocity dispersion 
and axis ratio respectively of th e projected SIE profile. As 
shown by lKormann et al.l (1 19941) . the lensing cross section for 
an SIE scales trivially with the Einstein radius Z?^ 



DlDls 



b = 4-7r 



(17) 



of the profile. Consequently, the distribution of halo orienta- 
tions for a lens sample, p(n) = CTB(n)/ {(Tb), is independent of 
the velocity dispersion ct,. of the halo. 

There is one last important element of the model that needs 
to be specified, namely the luminosity function of the sources 
being lensed. Here, we take the luminosity function to be a 
power law with slope of -2, which, while not exactly correct, 
is reasonably cl ose to the slope of the luminosity function of 
CLASS lenses (IChaell2003l: iMcKean et aTll2007l) . Moreover, 
this choice is ideally suited for numerical work since in such 
a case the biased cross section is simply ct^ = (/x) cr, implying 
that the biased cross section can be easily computed through 
uniform Monte Carlo sampling of the image plane. Since one 
of our goals in this work is to provide a test case for more 
complicated numerical algorithms, we choose a = -2. 

Having fully specified our model, we can now easily com- 
pute the biased lensing cross section for halos of any shape 
as a function of line of sight. Briefly, we proceed as follows. 
First, we compute the biased lensing cross section for SIE 
profiles as a function of the projected axis ratio q for a grid 
of q values. These data points are then fit using a third or- 
der polynomial fit, which we find is accurate to < 1%. Using 
this simple fit for (7h{q), and the fact that we can analytically 
compute the Einstein radius and projected axis ratio for a tri- 
axial halo along any line of sight, we can readily compute the 
mean lensing cross section of a halo averaged over all lines of 
sight. For a detailed description of our calculations, we refer 
the reader to the Appendices. 

Before we end, however, it is important to remark here 
that, despite its simplicity, we expect our model is more 
than adequate to investigate the qualitative trends that we 
would expect to observe in the data, and for providing order 
of magnitude estimates of the impact of triaxiality. Specif- 
ically, elliptical isothermal profiles appear to be excellent 
approximations to the true matter dist r ibution in re al lens 
systems (see e.g. iGerhard et al.l 120011: iRusin & Ma| 



2001 



2006 



Rusin et al."2003VRusin & Kochanek''2005'; iTreuetaD 
Koopmans et al. 2006; Gavazzi et al. 2007), so the triaxial 
isothermal mass distribution considered here should provide 
a reasonably realistic model for order of magnitude estimates. 
Wh ile more sophisticated m odels are certainly possible (see 
e.g. lJiang & Kochanekll2007l) . it is our view that the simplic- 
ity of the isothermal model more than justifies our choice of 
profile for a first pass at the problem. 

3.2. The Distribution of Halo Orientations for Triaxial 
Isothermal Profiles 

'By trivially, we mean a ocb-. 



Before we look at the distribution of halo orientations, it 
is worth taking a minute to orient ourselves in the coordinate 
system we have chosen. Consider first Eq. [15] The distance 
from the center of the halo to the intercept of a constant den- 
sity contour is maximized for the y axis, and minimized for 
the z axis, while the x axis is intermediate between the two. If 
we then parameterize the line of sight using the circular coor- 
dinates 9 and (p where 9 is the angle with the z axis and (p is 
the projected angle with the x axis, then our coordinate system 
is such that it has the following properties. 

• The x,y, and z axis of our coordinate system correspond 
to the middle, long, and short axis of the halo respec- 
tively. 

• Projections along cos(6') = 1 are along the short axis of 
the halo. 

• Projections along cos(6') = 0, = are along the middle 
axis of the halo. 

• Projections along cos{9) = 0, (/) = 7r/2 are along the long 
axis of the halo. 

The nice thing about this particular choice of coordinates is 
that in the cos(6') - cf) plane, both the long and the middle axis 
are represented by a single point, whereas the short axis is 
represented by an entire line. As we shall see, projections 
along the middle and long axis maximize the lensing cross 
section of a halo for quad and double lenses respectively, so 
having that maximum be a single point in the space of lines of 
sight is a desirable quality of our chosen coordinate system. 

Figure [T] shows the ratio b{fi) /bo where bo is the Einstein 
radius of an SIS with velocity dispersion cr,,, as well as the 
projected axis ratio q{n), for an isothermal ellipsoid with axis 
ratios qi = 0.75, q2 = 0.5. We can see the Einstein radius 
of the projected profile is maximized when projecting along 
the long axis of the halo, whereas the ellipticity is maximized 
when projecting along the middle axis of the halo, as it should 
be. Note we have only considered the range 9 E [0, 7r/2] and 
(pE [0, 7r/2] rather than the full range of possible lines of sight 
9 e [0,7r/2] and cp G [0,27r]. This is due to the symmetry of 
our model; all eight of the octants defined by the symmetry 
planes of the ellipsoids are identical. 

Let us now go back and study the distribution of line of 
sights for both doubles and quads. Figure |2] shows these dis- 
tributions for three types of halos: a prolate halo, an oblate 
halo, and a halo that is neither strongly oblate nor strongly 
prolate. As is customary, we parameterize the halo shape in 
terms of the shape parameter T which is defined as 



T = 



l-qj 
l-q2 



(18) 



Note that a perfectly prolate halo (^i = ^2) has T = I, whereas 
a perfectly oblate halo [qi = 1) has T = 0. From top to bot- 
tom, the halo shape parameters used to produce Figure |2] are 
T = 0.9 (cigar shape), T = 0.5 (neither strongly oblate nor 
strongly prolate), and T = 0.1 (pancake shape). The axis ratio 
q2 was held fixed at q2 = 0.5. Finally, the left column is the 
distribution of lines of sight for double systems, whereas the 
right column is the distribution for quads. For ease of com- 
parison, the color scale has been kept fixed in all plots. 

Let us begin by looking in detail at the doubles column first. 
As is to be expected, the distribution of lines of sight is peaked 
for projections along the long axis of the lens, as this line of 
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Fig. 1. — Top panel: Einstein radius h/b g as a function of line of siglit 
for a triaxial isothermal profile (see Eq. I15t with axis ratios q\ = 0.75 and 
qi = 0.5. Here, is the Einstein radius obtained in the spherically symmetric 
case = qi = 1). Bottom panel: Projected axis ratio for the same halo. Our 
coordinate system is such that cos(9) = 1 corresponds to projections along the 
short axis of the halo, cos((9) = 0, </< = corresponds to projections along the 
middle axis, and cos(9) = 0, = 7r/2 are projections along the long axis of the 
lens. Note projections along the long axis of the lens (cos(6) = 0,4) = 7r/2) 
result in the largest Einstein radius, whereas projections along the middle axis 
(cos(9) = 0, (/> = 0) maximize the ellipticity of the projected density profile. 



sight maximizes the Einstein radius of the projected profile. 
Moreover, the distribution is very sharply peaked for cigar- 
like halos (top row), but is rather flat for pancake-like halos 
(bottom row). The reason that the distribution of lines of sight 
for pancake-like halos is so flat is simple: for an oblate halo, 
projecting along either the long or medium axis of the halo 
results in a large Einstein radius, but also a large ellipticity, 
so a large part of the multiply imaged region of the source 
plane actually corresponds to four image configurations, tak- 
ing away from the cross section for producing doubles. When 
projecting along the short axis of the lens, the Einstein radius 
is minimized, but the projected mass distribution is nearly 
spherical, so the majority of the multiply-imaged region pro- 
duces only doubles. 

The column corresponding to quads has much more inter- 
esting structure. First, note that the distribution of line of 
sights for quad lenses peaks for projections along the mid- 
dle axis of the lens rather than the long axis of the lens. As 
noted earlier, projections along the middle axis of the lens 
maximize the ellipticity of the projected profile, so relative 
to projections along the long axis of the lens, it is evident that 



the increase in ellipticity more than offsets the slightly smaller 
Einstein radii for the purposes of enhancing the lensing cross 
section for producing quad systems. It is also interesting to 
note that while the peak of the distribution is always clearly 
about the middle axis of the lens, the shape of the distribu- 
tion varies considerably in going from prolate halos to oblate 
halos. In particular, note that for prolate halos the peak about 
the middle axis is relatively narrow. What is more, projections 
along the short axis of the lens are more likely than projections 
along the long axis because the latter minimizes the elliptic- 
ity of the projected profile. For oblate halos, on the other 
hand, projections along the long axis of the lens are almost as 
likely as projections along the middle axis. This is simply be- 
cause for such halos, there is little difference in the ellipticity 
of the projected profile between projections along the middle 
and long axis of the halos. Consequently, both axes result in 
highly effective quad lenses. Note too that for pancake-like 
halos, projections along the short axis are strongly avoided, 
since this projection minimizes both the Einstein radius and 
the projected axis ratio of the lens. 

In short, then, prolate halos and oblate halos will have very 
different orientation distributions: for prolate halos, nearly 
all doubles will be due to projections along the long axis of 
the lens, while most quads will be due to projections along 
the middle axis of the lens, followed by projections along the 
short axis. For oblate halos, however, all halo orientations are 
almost equally likely in the case of doubly imaged systems, 
whereas quads strongly avoid projections along the short axis 
of the halo. 

The remainder of the paper will explore whether these re- 
sults have a significant impact on the statistical properties of 
the halo population. Specifically, we will first consider the 
ellipticity distribution of lensing galaxies compared to that of 
galaxies as a whole. We will then discuss how these results 
affect the predicted quad-to-double ratio, and finally, we will 
investigate whether lensing halos are expected to have a sig- 
nificantly biased projected substructure mass fraction. 

4. THE PROJECTED AXIS RATIOS OF LENSING HALOS 

As mentioned in the introduction, if one assumes that the el- 
lipticity of the light and that of the mass are monotonically re- 
lated, then one would naively expect that lensing galaxies that 
produce quads ought to be more elliptical than the average 
galaxy because the lensing cross section for quads increases 
with increasing ellipticity. Similarly, galaxies that produce 
doubles should tend to be more spherical. In this section, we 
discuss the impact of halo triaxiality on the distribution of axis 
ratios for double and quad lenses. 

Given a line of sight n, we can compute t he ax is ratio ^(n) 
of the projected mass distribution (see Eq. IB24b . Using the 
distribution of lines of sight /c(n), one can then easily compute 
the distribution of projected axis ratios q for a sample of lenses 
via 

p{q\quq2) = J ^p(h)6D(q{n\quq2)-q)- (19) 

Figure [3] shows the distribution of the projected axis ratio of 
both quad and double systems for the sample pancake-like 
(oblate, r = 0.1) and cigar-like (prolate, T = 0.9) halos from 
Figure |2] As is to be expected, the distribution for quad sys- 
tems is considerably skewed towards high ellipticity systems, 
whereas the distribution for doubles is much flatter. More- 
over, the quads distribution is significantly more skewed for 
prolate (cigar-like) systems than for oblate (pancake-like) ha- 
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Fig. 2. — Contours of the oiientation distribution p(n) = erg/ (erg) for triaxial isotliermal profiles. The left and right columns show the distributions for double 
and quads respectively, while the three rows are for three different halos: the top row is for prolate (cigar-like) halos (T = 0.9), the middle row if a for a triaxial 
halo that is neither strongly oblate nor prolate (T = 0.5), and the bottom row is for oblate (pancake-like) halos (T = 0.1). The short to long axis ratio qi was 
held fixed at = 0.5. For ease of comparison, the color scale is the same in all plots. Note that for doubles, the distribution of lines of sight is peaked about 
projections along the long axis of the lens, as we would expect. For quad systems, however, the distiibution peaks for projections along the middle axis of the 
lens, which corresponds to maximizing the ellipticity of the resulting projected profile (see Figure[T). 
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Fig. 3. — Distribution of projected axis ratios for quads and doubles for a 
sample oblate halo (T = 0.1) and a sample prolate halo (T = 0.9). The distri- 
bution of projected axis ratios for doubles and quads are fairly different, with 
that of quads being peaked at low axis ratios (high ellipticities) while that of 
doubles is relatively flat. The quads distribution for prolate halos is signifi- 
cantly more peaked than that of oblate halos, suggesting that the distribution 
of projected axis ratios in quadruply imaged systems may help constrain the 
shape of the three dimensional matter distribution of halos. Also shown in the 
figure as lines on the top axis are the values of qojs, the projected axis ratio 
for which 75% of the lenses have q < cjojs (the oblate and prolate values for 
qoji are nearly identical). Note that the difference in go. 75 between quads and 
doubles is quite pronounced, with Aqojj > 0. 1 for both oblate and prolate 
halos. 

los. Based on Figure |3] we have attempted to distill the dif- 
ference between quads and doubles into a single number. We 
define the axis ratio qojs as the axis ratio for which 75% of the 
lenses have axis ratios q < (70.75 The value ^0.75 for quads 
and doubles for both sample halos is also shown in Figure |3] 
as lines along the top axis of the plot. It is clear that the pro- 
jected axis ratio ^0.75 for doubles and quads is very different, 
with A^o.75 > 0.1 for both oblate and prolate halos. 

Figure |4] shows the difference Aqojs between doubles and 
quads (i.e. qo'-jf''^" - qo'"s \ solid line) and between doubles 
and the overall halo populations (i.e. q'o'js^'^^ - q^'^'s ^ dotted 
line) as a function of the axis ratios q\ and q2. However, rather 
than using qi as an axis, we follow standard practice and pa- 
rameterize the shape of the halo in terms of the shape param- 
eter T defined in Eq. [18] There are several interesting things 
to be gathered from Figure |4] First, when comparing doubles 
to quads, note that while Aqoj<, is indeed large (^0.75 si 0.1) 
for both prolate and oblate halos, the difference can be larger 
for oblate halos than for prolate halos. Moreover, note that 
in going from oblate to prolate halos, the difference A^o.75 
goes through a minimum when q2 ~ q^ (solid line), in which 
case values as low as A^o 75 w 0.05 for ^2 ~ 0.5 are possible. 
Turning now to the comparison between doubles and random 
halos, we see that the difference in g'0.75 for these two halo 
populations becomes negligible in the case of oblate halos, 
reflecting the near uniform distribution of lines of sights for 
doubles for oblate halos (see Figure |2]i. On the other hand, 
the fact that most prolate doubles are seen along the long axis 
of the halo implies that Ago 75 between doubles and random 
halos must be significant, and thus doubles tend to be more 

The 75% number is selected in a somewhat ad hoc manner. Basically, 
we wanted qx to fall past the large prominent peak seen in Figure [3] and in 
that sense X = 80% or X = 90% would work just as well. On the other hand, 
observational estimates of qx for X close to unity would be quite difficult, so 
to some extent we wanted X to be as small as possible. We chose X = 75% as 
a reasonable value. 



Fig. 4. — Difference Aijo.75 between doubles and quads (q' 



,doiib!es 
0.75 



quads 
"^0.75 



solid line) and between doubles and the overall halo population (<?o 75*'''' — 
ijjj"'™, dotted lines) as a function of the small-to-large axis ratio qi and the 
shape parameter T (see Figure[3]for the definition of qojs)- For low values 
of q2 {q2 ~ 0.5), Ago. 75 between doubles and quads is typically large, of 
order 0.1. The minimum in Ago. 75 between quads and doubles occurs for 
axis ratios q2 ~ gj, the latter relation being shown above with the thick, solid 
straight hne. By comparing the axis ratio distribution for quad lenses to those 
of doubles and those of the overall halo population, one could in principle 
determine if halos are typically oblate (pancake-hke), prolate (cigar-like), or 
in between. 

circular than the typical halo. 

In short, then, the quantity Ago. 75 between doubles and 
quads and between doubles and random halos can, at least 
in principle, help determine whether most halos are oblate or 
prolate. If halos are prolate, the difference Ago. 75 between 
doubles and random halos is large. If this difference is small, 
we can then look at the difference Ago.75 between doubles 
and quads. If this last difference is large, then halos are typi- 
cally oblate, whereas if the difference is small, then halos are 
neither strongly oblate nor strongly prolate and q2^ q\- 

In practice, however, the above test is difficult to execute. 
In particular, while lens modeling can provide some mea- 
sure of the axis ratio q in quad systems, there remains a 
fair amount of uncertainty due to the approximate degener- 
acy between galax y ellipticity and external shear (see e.g. 
iKeeton et alj 119971) . This degeneracy is even stronger for 
doubly-imaged systems, and worse, there is no way of de- 
termining the axis ratio of the mass for non-lensing galaxies. 
Fortunately, at the scales relevant for strong lensing (< 5 kpc), 
bar yons dominate the total matter budget in early type galax- 
ies (iRusin et al.l l2003h . so one expects that the dark matter 
distribution in these systems will have the sa me ellipticity 
and orientation as the baryons. Observationallv. lKeeton et all 
(1998) (see also Keeton et al. 1997) compared the projected 
ellipticity of the light in lensing galaxies to the ellipticity re- 
covered from explicit lens modeling, and found that the light 
and the mass tend to be very closely aligned, though the mag- 
nitude of the ellipticities is not clearly correlated and the mod- 
est quality of the photometry available at the time made their 
elli pticity measuremen ts difficult. Moreover, the galaxy sam- 
ple lKeeton et al.l (Il998l) included many galaxies that had non- 
negligible environments that were not incorporated into the 
model. More recently, a detailed study of the Sloan Lens ACS 
Survey (SLACS Bolton et al. 2006) with more isolated galax- 
ies supports the hypothesis that the ellipticity of the light is in 
fact extremely well matched to the ellipticity of the projected 
mass, at least on scales comparable to the Einstein radii of 
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the galaxies dKoopmans et al.ll2006l) . ' ' Thus, for the purposes 
of this work, we simply take the isophotal axis ratio of lens- 
ing galaxies to be identical to the total matter axis ratio for the 
purposes of investigating whether lens biasing can be detected 
in current lensing samples. 

Figure|5]shows the cumulative distribution of isophotal axis 
ratios for quad lenses (solid) and double lenses (dashed) for all 
lensing galaxies in the CASTLES'" database with isophotal 
axis ratios measurements.'-' Of course, the selection function 
for this sample is impossible to quantify objectively, but our 
intent is simply to see whether any differences between lens- 
ing galaxies and random galaxies can be found. Also shown 
in the figure are the axis ratio distributions of early type galax- 
ies as reported b y two different groups: th e dotted line shown 
is the fit used bv lRusin & Tegmarkl(l2001l ) to model the distri- 
bution of axis ratios in early type ga laxies based on measure- 
ments by iJorgensen & Franxl (1 19941). and is also quit e close 
to the distribution recovered bv iLambas etaP (11991) . The 
dashed-dotted li ne is the axis ratio distribution obtained by 
iHao et al.1 (l2006l) using the SDSS Data Release 4 photometric 
ca talog, and is a very clo s e match to the dis tribution recovered 
bv lFasano & \^ (Il991l) . iHao et al.l (l2006l) note d that it is un- 
clear w hy these two distributions differ, though iKeeton et aP 
(|1997|) note that such a difference can easily arise depending 
on whether SO galaxies are included in the galaxy sample or 
not (with SO galaxies being more elliptical). Here, we simply 
consider both distributions. 

Given that the axis ratio distribution for both quads and 
doubles largely fall in between the two model distributions 
we considered, it is immediately obvious that no robust re- 
sults can be obtained at this time. Specifically, uncertainties 
in the details of the selection function of the galaxies used to 
construct the isophotal axis ratios are a significant systematic. 
More formally, using a KS-test, we find that the isophotal axis 
ratio distributions of both quad and double lens galaxies are 
consistent with that of the early type galaxy population as a 
whole (irrespective of which model distribution we choose) 
and with each other as well. Interestingly, whether or not 
we restrict ourselves to galaxies that are isolated or whether 
we include all lensing galaxies does not appear to change 
the result in any way. Naively, then, the consistency of the 
axis ratio distributions suggests that halos are typically neither 
strongly oblate nor prolate, but rather somewhere in between, 
where the quantity A^o.vs exhibits a minimum, which occurs 
at qi ~ q\. 

Given that current lens samples are too small for detecting 
any difference on the ellipticities of quadruply and doubly im- 
aged systems, it is worth asking whether or not a detection is 
possible in principle. That is, how many lenses must one have 
in order to detect quad systems as being more elliptical than 
doubles? To answer this question, we need to first assume a 
simple model for the distribution of axis ratios Q{q\, qi), with 
which one could then compute the resulting projected axis ra- 
tio distributions for doubles, quads, and the galaxy population 
at large. We should note here, however, that in detail our re- 
sults will depend on the adopted distribution Qiqi^qi), which 
is not known. 

It is not immediately obvious what the most correct model 

' ' We note, however, that SLAC lenses tend to have Einstein radii that are 
quite comparable to their optical radii, so the agreement is really expected. 
In principle, a discrepancy could exist for lenses with larger Einstein radii for 
which the total mass has a larger dark matter component. 
http://cfa-www.harvard.edu/castles/ 

This data was kindly provided by Emilio Falco, private communication. 
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Fig. 5. — Cumulative distribution of the isophotal axis ratio for early type 
galaxies. The step-lines shown correspond to galaxies that produce quads 
(solid) and doubles (dashed). The smooth curves show the overall early type 
galaxy population in the Coma clust ers (dotte d, see Rusin & Tegmark 200ll) 
and in the SDSS DR4 (dashed, see IHao eTa l. 2006). The axis ratio distri- 
bution for lensing galaxies is consistent with both model distributions, and 
with each other, suggesting that the typical halo is neither strongly oblate nor 
prolate, with axis ratios qi and q2 such that qi^ q\. 

distribution Q{qi,q2) should be. While there have been many 
studies that have investigated the distribution of axis ratios 
of dark matter halos in simulations (see e.g. [ Warren et al.l 
fT99llJing & Sutoll2002HBailin & Steinmetzll2005h . it has be- 
come clear that the distribu tion itself depends on many vari- 
ables , including halo mass dKasun & Evrardll2005l: iBett et all 
12007b, radius at wh ich the shape of the halo is measure d 
dHavashi et"aLl l2007h , halo environment (iHahn et al.l l2007l) , 
and whether the halo under consideration is a p arent halo or 
a subhalo of a larger object (iKuhlen et al.ll2007l) . Adding to 
these difficulties is the fact that different authors use different 
definitions and methods for measuring the shapes of halos, 
which forces one to go to great lengths in order to ensure a 
fair com parison of the result s from different groups (see for 
example lAllgood et al.ll2006l) . Even more problematic that all 
of these difficulties, however, is the fact that not only can the 
distrib utions of baryons have a d ifferent shape from the dark 
matter ( iGottlober & Yepesll2007l) , baryons dominate the mass 
budget in the halo regions where strong lensing occurs, and 
can therefore drarnatical l y impact halo shapes at those scales 
(iKazantzidis et al.1 120041: iBailin et alll2005l: iGustafsson et all 
|2006^ . Since our intent here is simply to provide a rough 
estimate of the number of lenses required to detect a signif- 
icant difference in the ellipticities of quad and double sys- 
tems, we simply adopt a fiducial mode l that is based primarily 
on the results of I All good et al.l ( |2006|) and IKazantzidis et al.l 
(2004), and use it to estimate the number of lenses neces- 
sary to detect the lar g er elli pticity of quad systems. Specif- 
ically, lAUgood et"an (|2006|) obtain that for an M, halo the 
distribution of the short-to-long axis ratio of dark matter ha- 
los is Gaussian with a mean of ( q^) = 0.54 and a s t andar d 
deviation a^^ =0.1. As noted bv IKazantzidis et akl d2004l) , 
baryonic cooling tends to circularize the mass profiles of ha- 
los, so we adopt instead a somewhat larger ratio ((72) = 0.65, 
but retain the dispersion a^, =0.1. The adopted value for 
(q'2) is larger than that obtained from dissipationless sim- 
ulations, but sma l ler th an that found in the simulations of 
IKazantzidis et al.l d2004l) . as the latter suffer from the well 
known over-cooling problem and therefore overestimate the 
impact of baryons on the profiles. In addition, we truncate 
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Fig. 6. — Cumulative distribution of the isophotal projected axis ratio for all 
galaxies, quad lenses, and double lenses, using the fiducial model discussed in 
the text. For comparison, we also show the o bserved distribution of isophotal 
axis ratios in the SDSS DR4 as measured by Ha(3et al. ('2()6 6[). Keep i n min d 
that the model is not meant to be a fit to the data from H ao et al] 120061) , 
but rather is only meant as a useful approximate model that allows us to 
estimate the number of lenses necessary to distinguish between the ellipticity 
distribution of quads and doubles. We find that N 350 lenses are needed to 
empirically distinguish the two distributions at the 95% confidence level. 

the distribution at q2 = 0.4, as the expressions for the lensing 
cross sections are no longer valid for systems with projected 
axis ratios below 0.4,''* Finally, the val ue of the intermediate 
axis qi is obtain following th e model of lAllgood et al.l (|2006|) 
(itself based on the work bv l.Ting&Sutoll2002h . namely, the 
quantity p = q2/q\ is drawn from the distribution 



p{p\s) = 



2(1- s) 



1- 



2p-\- 



(20) 



where s = min(0,55,^2)- 

Figure|6]shows the cumulative distributions of the predicted 
isophotal axis ratios for all galaxies, as well as for quad and 
double systems. Also shown for referenc e are the axis ra- 
tio measurements of early type galaxies by iHao et al] (l2006h 
using SDSS DR4 data. Note that, as we expected, the dif- 
ference in the axis ratio qojs between doubles and quads is 
of order 0.05, The maximum vertical distance between the 
cumulative distributions functions for quads and doubles is 
D K, 0,15, which, using a KS-statistic, implies that roughly 
300 lenses (150 quads, 150 doubles) with good isophotal mea- 
surements are necessary to detect the difference between the 
two distributions at the 95% confidence level, A 5cr detection 
would require « 1,400 lenses. Such large number of lenses 
is larger than the current list of known lensing systems, but is 
certainly within the rea lm of what one rnay expect from future 
lens searches (see e.g, iKoopmans et alJl2b04t [Marshall et alJ 
l2005h . 

5, TRIAXIALITY AND PREDICTIONS FOR THE QUAD-TO-DOUBLE 

RATIO 

We showed above that halo triaxiality can have an important 
impact on the distribution of axis ratios for lensing galaxies. 

For SIE profiles, if the projected axis ratio q < 0.4, then naked cusp 
configuration appear. Since the analytical formulae we used to compute a{q) 
all compute the area contained within the tangential caustic, it follows that 
for q < 0.4, our cross section estimates would correspond to the total cross 
section for producing either quads or naked cusps. To avoid this complica- 
tion, we simply truncate our axis ratio distribution at q2 = 0.4. Note however 
that since qi = 0.4 is already 2.5cr away from the adopted mean we expect 
the introduced cutoff to have a negligible impact on our results. 



Since the projected axis ratio of a halo plays a key role in the 
expected quad-to-double ratio of lensing galaxies, it is easy 
to see that triaxiality should also affect this statistical observ- 
able. This is the problem we wish to consider now: how does 
triaxiality affect the quad-to-double ratio of lensing galaxies? 

Consider first equation [TT] For our semi-analytic case, the 
halo parameters P that determine the mass distribution of the 
halo are simply the halo velocity dispersion and its two axis 
ratios q^ and q2- What is more, we saw that if we define i'o(c,,) 
as the Einstein radius of an SIS of velocity dispersion (jy, then 
the ratio (j/b^ depends only on the axis ratios qi and q2. If 
we make the further assumption that the distribution of halo 
parameters is separable, i,e, that 



drihah 



dill, 



lalos 



dzhd(Jydq\dq2 dzhdcTy 



Q(qi,q2), 



(21) 



then it is easy to see that the ratio of the total number of quad 
systems to double systems depends only on the distribution 
of axis ratios Q{qi,q2) because the overall scaling of the lens- 
ing cross sections for both doubles and quads just factors out 
of the problem. Thus, the ratio of quad-to-doubles is given 
simply by 



riquqi)-- 



No, of quads ys mAi 



,(4) I 



No, of doubles 



(2)1 



(22) 



The top panel of Figure|2]shows the dimensionless mean bi- 
ased lensing cross section (ctb) /^o for both doubles and quads 
averaged over lines of sight for a population of randomly ori- 
ented halos. Also shown in the bottom panel is the quad- 
to-double ratio. As expected, large (> 0,3) quad-to-doubles 
ratios require strong deviations from spherical symmetry, so 
q2 needs to be small. Interestingly, however, all of the con- 
tours in both the top and bottom panel of Figure |7] are nearly 
vertical: lensing cross sections are nearly independent of halo 
shape. We can understand this qualitatively as follows. In the 
case of doubles, there is a tradeoff between two competing 
effects: for 1 > Q'l ^ q2, there are many lines of sight that en- 
hance the Einstein radius of the lens, but only moderately so. 
For \':$> q\ > <?2 on the other hand, there are only a few lines 
of sight that enhance the Einstein radius of the lens (i,e, pro- 
jections along the long axis of the halo), but the enhancement 
is much greater. Thus, the overall boost to the Einstein radius 
is offset by the reduced "volume" of lines of sight available 
for forming doubles and vice versa, A similar effect occurs 
for quads: oblate halos make effective lenses when projected 
along either the long or middle axis of the lens, but strongly 
avoid the short axis, so the "volume" of lines of sight available 
to oblate halos is small. Prolate halos, on the other hand, are 
not quite as effective as oblate halos at making quads, but can 
produce quads over a larger range of possible lines of sight. 

At any rate, one thing that is clear from Figure |7] is that 
halo shape does not have a significant impact on the expected 
quad-to-double ratio. One extremely interesting consequence 
of this results is that it implies that halo triaxiality can be 
properly incorporated into lensing statistics studies without 
greatly increasing the number of degrees of freedom in the 
problem. More explicitly, traditional lens statistics studies 
use as input the observed two dimensional ellipticity distri- 
bution of early type galaxies, and approximate the effects of 
triaxiality by multiplying the usual isothermal ellipsoidal pro- 
files with a normalization factor computed assuming hal os are 
either all perfectly oblate, or perfectly prolate (see e,g, IChad 



10 



ROZO ET AL 




0.6 O.i 
q2 [Short/Long] 



1.0 




0.6 0.7 O.i 
q2 [Short/Long] 

Fig. 7. — Top panel: The dimensionless biased lensing cross section (Je/b^ 
for quads (solid) and doubles (dotted) as a function of the axis ratios q[ and 
q2 of the halo population. Bottom panel: quad-to-double ratio as a function 
of the axis ratios q\ and q2 for triaxial isothermal profiles. We only consider 
the range q2 > 0.4 because below this value our semi-analytical calculations 
based on Kormann et al. ( 1994) break down. Note the observed ratio of about 
0.4 can only be obtained for halos that deviate strongly from spherical sym- 
metry. Interestingly, all contours in both panels are very nearly horizontal, so 
whether halos are prolate (cigar-like) or oblate (pancake-like) has almost no 
impact on the lensing optical depths or the quad-to-double ratio. 



12001 120071: IOgura2007h . The main reason this is done, rather 
than considering triaxial halos and averaging over lines of 
sight, is that in order to do the latter calculation, one needs 
to know something about the distribution of axis ratios. We 
have shown, however, that such a calculation would in fact be 
nearly independent of assumptions made about the intermedi- 
ate axis qi. In other words, a proper calculation that weights 
lines of sight according to their biased lensing cross section 
rather than uniform weighting (as implicitly done when tak- 
ing the ellipticity distribution to be that of early type galax- 
ies as a whole) effectively involves no more freedom than the 
usual approach, the main difference being that the assump- 
tions made will involve not the ellipticity distribution, but 
rather the distribution of the short-to-long axis ratio q2, which 
can itself b e constrained using the projected ellipticity distri- 
bution (e.g. lLambas et al.|[T992l) . 

6. THE SUBSTRUCTURE MASS FRACTION IN THE INNER 
REGIONS OF LENSE^G HALOS 

One of the important predictions of the CDM paradigm 
of structure formation is that galactic halos contain a 
large amount of bound substructure within them (see e.g. 
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Fig. 8. — The projected substructure mass fraction within 3% of the virial 
radius as a function of cos(6;o„g ) where Biong is the angle between the projec- 
tion axis and the long axis of the halo. En'or bars show the dispersion in the 
sample rather than the error on the mean. The solid line shows a simple con- 
stant plus gaussian fit of the form /,(.v) = 0.003 + 0.018exp(-(l -.vj'/O.l^) 
where x = cos((?/o„g). Note projections along the long axis of the lens have 
significantly higher substructure mass fractions than the typical projection. 



IWhite & Reeslll978l: iBhimenthal et allll984l) . Observation- 
ally, however, both our own galaxy and M31 have an or- 
der of magnitude less luminous companions than is predicted 
if o ne assumes substructur e s have a fixed mass to light ra- 
tio (iKauffmann et al.|[T99l iKlvpin et all 1 19991: iMoore et al.l 
119991) . Currentlv. the favored explanation for this discrepancy 
is that the mass to light ratio of such small structures depends 
strongly on the history of the objects, and therefore only a 
select subse t of the su bstructures within the halo become lu- 
minous (e.g.lSomervill e & Primack 1999l: lBenson et al.l2002l : 
iKravtsov et al.ll2004tlSales et al.ll2007ft . While such scenarios 
appear to be in good agreement with the data, it would still 
be desirable to provide as direct detection as possible of the 
remaining dark substructures. 

Motivated by the fact that dark substr uctures can only be 
discovered via their gravitational signal, iDalal & Kochanekl 
( t2002 al) investigated whether the well known flux anomalies 
problem could be explained as the action of dark substructures 
embedded within the halo of the lensing galaxy. Using a rel- 
atively simple model, they found that in order to explain the 
observed flux anomalies, one requires a projected substruc- 
ture mass fraction in the range 7% > fs > 0.6% at the 90% 
confidence level. It was then argued bv lMao et al.l ( |2004|) that 
such a substructure mass fraction was slightly larger than the 
mass fraction obtained from simulations /, « 0.5%. 

Recently, it has become clear that the distribution of sub- 
structures in dark matter halos is not spherically symmetric, 
but is instead triaxial, and aligned with the major axis of the 
halo. Since lensing halos are not randomly oriented in space, 
the mean projected substructure mass fract i on for all halos - 
the f, w 0.5% value obtained bv lMao et all (12004 - need not 
be the same as the mean substructure mass fraction for lens- 
ing halos, which would in turn affect the oretical predictions 
(e.g. iRozo et al.|[2006bl: IChen et al.|[2007l) . Here, we use the 
results on substructure alignments in numerical simulations 
to estimate the dependence of the projected substructure mass 
fraction /, on the projection axis. More specifically, assuming 
that substructures do not significantly alter the biased lensing 
cross sections for the halos, we compute the mean substruc- 
ture mass fraction for doubles and quad lenses as a function 
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Fig. 9. — Mean projected substructure mass fraction (fslqi^cji) as a func- 
tion of halo axis ratios for doubles (dotted lines) and quads (solid lines), 
shown above as percentages (i.e. we show 100 ,52))- The projected 
substructure mass fraction as a function of line of sight is taken to be the fit 
shown in Figure[8] and biased cross section weighting assumes substructures 
do not significantly impact the lensing cross section. For comparison, the 
mean substructure mass fraction for randomly oriented halos in our model is 
(/,) = 0.46%. 

of the axis ratios qi and q2 of the lensing halos. 

We begin by presenting the substructure mass fraction 
as a function of line-of-sight in simulated dark matter ha- 
lo s. In Figure [8] we reproduce the distribution as presented 
in IZentnen (l2006h . This figure shows the mass fraction pro- 
jected within 3% of the virial radius as a function of the pro- 
jection angle cos{6io„g) for a sample of halos in a dissipa- 
tionless A^-body simulation of structure growth. The angle 
Oiong is defined as the relative angle between the projection 
axis and the long axis of the halo. The data for the figure 
come from 26 host dark matter halos with masses in the range 
IO'^/i-'Mq <M< lO'^/r'M0, and the error bar shown rep- 
resents the dispersion in the sample rather than the error on 
the mean. The halos were drawn from a high-resolution flat, 
ACDM simulation with rj,„ = 0.3, ag = 0.9, ^bh^ = 0. 023, and 
h = 0.7. Details on the simulations can be found in IZentneJ 
(I2OO6) or in Gottlober & Turchaninov ( 2006). 

Using the fit to fs{cos{6io,jg)) shown in Figure [8] we com- 
pute the mean projected substructure mass fraction (/, \qi , ^2) 
for a population of double and quad lenses as a function of 
the halo axis ratios qi and q2- Our results are shown in Figure 
|9] For reference, the mean substructure mass fraction for ran- 
domly oriented halos obtained from the fit shown in figure |8] 
is (/s) = 0.46%. As per our expectations, we find that prolate 
(cigar-like) doubles have substructure mass fractions that are 
enhanced relative to the average halo, with (/j) w 0.6%. Note 
though that this enhancement is relatively minor, and slowly 
decreases to the random average as halos become oblate. 

More interesting to us is the behavior of quads, for which 
we find a mild enhancement relative to random for oblate ha- 
los, and a decrease in the expected substructure mass fraction 
for prolate halos. This can be easily understood from Figure|2] 
oblate quads strongly avoid projections along the short axis of 
the halo, and projections along the middle and long axis of the 
lens are nearly equally likely. Consequently, one expects an 
enhancement of the substructure mass fraction because some 
lines of sight with low /, are avoided. On the other hand, for 
prolate halos, projections along the long axis of the lens are 
the least common, so indeed we expect the mean projected 
substructure mass fraction for these systems to be reduced. 



Overall, though, it is clear that for quad systems - which 
are the only kind of sy stems for which ( f, ) may b e estimated 
using the methods of iDalal & Kochanela (l2002ah - the sub- 
structure mass fraction in the inner regions of a halo can- 
not be significantly enhanced due to lens biasing if the im- 
pact of substructures on the lensing cross section of galac- 
tic halos can be neglected. Thus, lens biasing does little to 
soften the slight (and in these authors' opinion, not terribly 
si gnificant) discrepancy betw een the values of recovered 
bv lDalal & Kochanekl ( l2002ah and those from numerical sim- 
ulations. 

7. CAVEATS AND SYSTEMATICS 

Before we finish, we believe it is important to mention two 
systematics that could significantly affect the conclusions pre- 
sented in this work. Specifically, throughout we have assumed 
that the lensing cross section is dominated by the smooth mass 
distribution of lensing galaxies, and we presented in section[3] 
several studies that suggest that our model for the mass dis- 
tribution of early type galaxies is a reasonable one. As men- 
tioned in the introduction, the possible discrepancy between 
theory and observation concerning the quad-to-double ratio of 
the CLASS lenses has raised the possibility that lensing cross 
sections are in fact heavily influenced by the environment of 
the halo or possibly by substructures with in it. We briefly 
discuss each of these in turn. 

We begin by discussing halo environments. In our calcula- 
tions above, and in most of the lensing statistics literature, the 
effect of halo environment on lensing statistics is neglected. 
This is not an entirely ad hoc assumption. Theoretical es- 
timates of the amount of shear that the typical lens experi- 
ences are quite small (7 « 0.02, see e.g. iKeeton et alj|1997l : 
iDalal & WatsonI l2004t) . so its impact should be negligible. 
Curiously, however, explicit lens modeling of known systems 
usually requires large external shears (7 « 0.1 ) in order to 
provi de reasonable fits to observations (see e.g. iKeeton et alj 
Il997h . Moreover, direct estimates of the environment of lens- 
ing galaxies also s upport the idea o f a stronger effect from 
nearby structures ( lOguri et al.l l2005h . The discrepancy be- 
tween these observations and the predictions for halo envi- 
ronments are themselves an interesting problem, which ul- 
timately may or may not be related to the usual quad-to- 
double ratio problem. At any rate, one might hope that 
even if such large external shears are correct, their impact 
on the quad-to-double ratio would still be negligible if their 
orientation is random. This expec tation was indeed con- 
firmed by ( iRusin & Tegmarlj 1200 Ih . Unfortunately, it is 
known that a significant fraction of lenses a re actually mem- 
ber galaxies of intermedia te mass groups dMomcheva et alj 
I2OO6I: IWilhams et alj|2006l) . and that galaxies in groups and 
clusters tend to be radially aligned (Pereira & Kuhn 2005!; 
iDonoso et aP l2006t iFaltenbacher et alJ l200% . implying the 
randomly oriented shear assumption is Ukely not justified. In- 
deed, a careful analysis of the impact of the halo environment 
for group members shows neglecting to take said environ- 
ment into account can lead to an underestimate of the ratio 
of the quad-to-double lens ing cross sections for such galaxies 
as large as a factor of two (IKeeton & Zabludoffl2004) . At this 
point, what seems clear is that there is not as of yet a defini- 
tive answer as to exactly how important galaxy environments 
are, and thus, we have opted for making the simplest possi- 
ble assumption for the purposes of this work, that is, we have 
ignored the impact of large-scale environments. 

The second solution to the quad-to-double ratio problem in- 
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volves substructures. Specifically. ICohn & Kochanekl ( |2004|) 
have shown that the lensing cross section of galaxies is 
severely affected by substructures. If this is indeed the case, 
the way in which lensing galaxies are biased relative to the 
overall galaxy population depend not only on the smooth 
component of its mass distribution, but also on the spatial 
distribution of substructures within the galaxy halo. Interest- 
ingly, in such a scenario halo triaxiality would impact the ori- 
entation of halos relative to the line of sight now only through 
the biasing due to the smooth matter component, but also be- 
cause of the previously mentioned alignment between the sub- 
structure distributions and the smooth mass distributions. We 
leave the question of exactly how such a population of halos 
would be biased to future work (Chen et al., in preparation). 

8. SUMMARY AND CONCLUSIONS 

The triaxial distribution of mass in galactic halos implies 
that the probability that a galaxy becomes a lens is depen- 
dent on the relative orientation of the galaxy's major axis to 
the line of sight. Consequently, a subsample of randomly ori- 
ented galaxies that act as strong lenses will not be randomly 
oriented in space. The relative orientation and the strength of 
the alignment depends on the shape of the matter distribution, 
and on the type of lens under consideration: prolate doubles 
have a high probability of being project along their long axis, 
whereas the distribution of oblate doubles is nearly isotropic. 
Prolate quads are most often projected along their middle axis, 
though the degree to which alignment occurs is not as strong 
as for prolate doubles. Interestingly, highly prolate quads are 
also more likely to be projected along their short axis than 
along the long axis, though this very quickly changes as halos 
become more triaxial and less prolate. Oblate quads strongly 
avoid projections along the short axis of the lens, but projec- 
tions along the other two axis are almost equally likely. 

An important consequence of the differences in the distri- 
bution of halo orientations for quad lenses, double lenses, and 
the galaxy population as a whole is that the ellipticity distri- 
bution of these various samples must be different, even if the 
distribution of halo shapes is the same. Specifically, we pre- 
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diet that quad lenses are typically more elliptical than random 
galaxies, and that the ellipticity distribution of doubles is very 
slightly more circular than that of random galaxies. While 
current data do not show any indication of these trends, we 
have shown that w 300 (1 ,400) lenses are necessary to obtain 
a 2(7 (5(t) detection of the effect. 

The fact that halo triaxiality affects the ellipticity distribu- 
tion of lensing galaxies also means that halo triaxiality needs 
to be properly taken into account in lensing statistics. Conse- 
quently, we estimate how the biased lensing cross sections of 
galaxies depend on halo shape, and find that they are nearly 
independent of the halo shape parameter T. Instead, the mean 
biased cross section of a lens depends almost exclusive on the 
distribution on the short-to-long axis ratio (often denoted 
by s). 

Finally, given that the distribution of substructures in nu- 
merical simulations is observed to be preferentially aligned 
with the long axis of the host halos, we estimate how the pre- 
ferred orientation of lensing galaxies affects their predicted 
substructure mass fraction. We find that biases due to non- 
isotropic distribution of halos relative to the line of sight have 
an insignificant impact on the mean substructure mass frac- 
tion of lensing galaxies. 
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APPENDIX 

LENSING CROSS SECTIONS OF SINGULAR ISOTHERMAL ELLIPSOIDS 

The Singular Isothe rmal Ellipsoid (SIE) is one of the simplest lens models that can produce quadruply imaged sources. 
iKormann et al.l (1 19941) performed a detailed study of the lensing properties of SIE lenses, and, in particular, derived simple 
expressions for the total area contained within the tangential and radial caustics of such lenses. Specifically, given an SIE profile 

Jqal 1 
20 jc^ -I- q'^y'^ 

iKormann et al.l (Il994h found that the area (7r and a, contained inside the radial and tangential caustics is given by 

4-q cos~'(x) 



and 

a,= —^ I dx I • ' " -cos"'(x) I ^ " , (A3) 




respectively. Moreover, they showed that for q > qc where q^ w 0.394, the tangential caustic is entirely contained within the 
radial caustic, and hence the lensing cross section ct***' for forming four image lenses is simply cr*"^* = a,. Likewise, the lensing 
cross section for forming doubles is given by cr*^^ = (T, -cr,. 

Unfortunately, as derived in section 12.11 the relevant quantity for lensing statistics of a flux limited sample is not the lensing 
cross section itself, but the biased cross section as- Moreover, the latter cross section requires one to compute the magnification 
distribution for double and quad lenses, for which there are no closed form expressions. In this appendix, we numerically 
compute the magnification distribution pifi), and its first moment (/x) for both doubles and quads, and use them to compute the 
biased lensing cross section ctb = {fi) a appropriate for a source luminosity function ns{L) cx L~^. 
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Fig. 10. — Left panel: Distribution of magnifications p(p) for double and quadruple lenses for an SIE profile with axis ratios q = 0.4 and q = 0.8. The q = 0.4 
distribution for doubles is shown as a solid hne, and the corresponding distribution for quads is shown as diamonds. The dotted curve and trian gles are the double 
and quad distiibutions for q = 0.8. The solid curves on top of the diamonds and tiiangles are our empirical fit for quad lenses (see Eq. IA4I and is accurate to 
fa 5%. Right panel: Biased lensing cross section = {/J,) o" as a function of axis ratio for SIE profiles and assuming a source luminosity function with a power 
law slope a = -2. The solid curves are simple cubic fits which are accurate to better than 1%. 

The left panel of Figure[lO]shows the magnification distribution for doubly and quadruply image systems for SIE profiles with 
axis ratios q = 0.4 and q = 0.8. Note that the magnification distribution for doubles is very rich in features. The magnification 
distribution for quads, on the other hand, is relatively simple, and we can provide a simple fitting formula for it. To do so, first 
note that we know that in the limit fi oo, p'-'^\^i) oc /z~ , so we expect that p'-^\fi) w Nx~^f{x) where x = fi/fimin and /i„„„ is 
the minimum magnification for quad lenses, is a normalization constant, and f(x) is a function which asymptotes to unity and 
deviates from unity only for x w 1 . Consequently, we expand f{x) in a power series in terms of x"' , of which we expect only the 
first few terms would be necessary to produce a good fit. As it turns out, we found that f{x) needs only one non-constant term to 
result in excellent fits to p^'^\pL), and our final fitting function for p^'^\fi) is thus 



1 



(4) 



\+a/2 



^(l+fljc--). 



(A4) 



A priori, we would expect that the best fit value of the a coefficient in the above expression would be a function of the axis ratio 
q of the profile. While there does appear to be some such dependence, it is extremely mild, so we have opted for keeping a fixed 
to the value a = 0.83. We found that this expression is accurate to better than 5% for /i,„/„ < /i < 2QiJ.,„i„ and q > 0.4. 

The right panel of Figure [TO] shows the actual quantities we are interested in, the biased lensing cross sections (Jb = (/i) cr. As 
is obvious from the figure, the form of these biased cross sections is very simple, so even a simple quadratic fit results in quite 
good fits (of order a few percent). Since we wish our empirical fit to be accurate, we fit the numerically computed cross sections 
with a cubic, which is enough to obtain sub-percent level accuracy. Our best fit curves (in a least square sense) are 



4^' = -6.902 + 42.937^-33.240^^9.736^^ 
20.833^+13.236^^-3.816^1 



cr*f*= 11.409- 



(A5) 
(A6) 



Of course, we could have just as easily splined the numerically estimated values to compute the lensing cross section at any 
axis ratio q. We opted to fit the cross sections with a simple form both for simplicity, and in the chance that the fitting formulae 
provided here will be useful for other works. 



PROJECTED SURFACE DENSITY PROFILES OF TRIAXIAL ISOTHERMAL HALOS 



Consider an SIS profile 



Its triaxial generalization takes the form 



Psisir) ■■ 



1 



2itG r- 



(Bl) 



PsiE{^)=N(qi,q2) 



2^Gx^/q\ + f + Z^/ql 



(B2) 



where qi and q2 are the halo's axis ratios, and we have chosen a coordinate system that is aligned the halo's principal axis and 
such that \>q\> qi- p2 remains the ratio of the small to large axis. The prefacto r N(qi , q2) represents a relative normalization 
for halos of varying axis ratios which we will compute shortly. First however, since iKormann et aH (1 19941) use the notation where 
the axis ratios multiply rather than divide the coordinates, we rewrite the mass density as 



PsiE(^) = N(quq2) 



1 



27rG p\x^ + p^y^+z'' 



(B3) 
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where pi = q2/q\, pi = qi, and N{pi,p2) = q\N{q\,q2)- Note pi is the ratio of the small to middle axis, while pi remains the 
ratio of the small to large axis. We choose the normalization function Nipi^pi) such that the mass contained within a radius 
r is independent of the axis ratios p[ and p2, as appropriate if one wishes to investigate the impact of triaxiality on lensing 
cross sections at fixed mass with the latter defined using spherical overdensities. Integrating the above profiles and setting 
Msis(r) = MsiE(r) results in'^ 



N{puP2) = 



7r/2 



tan 



\/a(\—a) 



(B4) 



where we have defined a{(f);qi,q2) via 



a((l);qu qi) = p\ cos^(0) + p\ sin^(0). 



(B5) 



We wish to project psiE along an arbitrary line of sight. Let x be a coordinate system such that the z axis is aligned with the 
line of sight. We choose the x and y axis to be such that a rotation by an angle 6 along the y axis followed by a rotation along the 
z axis by an angle recovers the coordinate system x from Eq. IB3I The corresponding rotation matrix is given by 



R = 




sin 6' cos (; 
cos 6' 



By construction, the corresponding projected surface density is given simply by 

/•OO 

S(jc,y)= / dzpsiEiR'^) 



which has the form 



where 



^(x,y) = N{puP2) 



2ttG 



dz 



1 



A+Bz+Cz- 



and 



C = p\ sin^ 9 cos^ 4'+P2 sin^ sin^ (/> + cos^ ^ 



2 2 2 2 2 2 2 

A.o- = /5[Cos 0COS 0+/72COS 6* sin (/) + sin 
A = sin(2(/)) cosC^X-p 1 + pi) 
Ayy = p\ sin^ 0+ p\ cos^ 
By^ = sin(20)(/7 1 cos^ </> + P2 sin' </> - 1 ) 



By = sin(9) sin(20)(-/7i + pj). 



Note that if ^1 = ^2 = 1, then A„ = Ayy = C = I and Ai,. = Bx = By = 0, exactly as it should. Performing the integral in Eq. 
find 

1 



(B6) 



(B7) 

(B8) 

(B9) 
(BIO) 

(BID 

(B12) 
(B13) 
(B14) 
(B15) 
(B16) 
we 



^{x,y)=N(pi,p2) 



which has the generic form 



where 



i:(x,y) = N(pi,P2):r^ 



2G AC -By 4 
1 



2G (a^xX^ + a„xy + ayyy'^y 



a,,=A„C-Bl/4 
a^=AxyC-B,By/2 

ayy=AyyC-By4. 



(B17) 



(B18) 



(B19) 
(B20) 
(B21) 



For q\= q2 = 1, the above expressions reduce to axx = ctyy = 1 and a.vy = as appropriate for an SIS profile. For the more general 
case it is evident from equationuation IbTS] that using an additional rotation of the x-y plane we can diagonalize the projected 
mass density Yi{x,y). We find that the required rotation angle ?/; is given by 



tan2V'= — — . 

" To obtain the expressions above, we perform first the radial integral and then the integral where 9 is the azimuthal angle. 



(B22) 
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Using a ~ to denote the new coordinate system, we can thus write 



where 



q^ = ^ (B24) 

axx 

CTy = — j=^^a,, (B25) 

axx = ctxx cos^ ip + a vj. sin ■0 cos -0 + a,.,, sin^ ip (B26) 

a^y = axx sin^ -tp - sin ip cos -tp + ayy cos^ ip. (B27) 

As expected, the above expression for q reduces to ^ = pi/pi = q\ when we project along the z axis (i.e. the short axis), to q = p2 
when projecting along the y axis (i.e. the long axis), and to q = pi= q2/q\ when projecting along the x axis (i.e. t he middle axis). 
The pa rticular form of the parameterization of the surface density in Eq. IB23l is meant to match the conventions in lKormann et alJ 
(1 19941) . which was chosen to ensure the mass contained within a given density contour be independent of q for fixed ct,,. 
and the 



and we have defined 



